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ABSTRACT 


Fast automatic mesh generation schemes are required in many 
PPlication areas such as finite element analysis. The present 
'ork outlines one such methodology of mesh generation that can be 
mplemented on a parallel architecture of processors. The proposed 
lesh generation algorithm uses a modified quadtree approach. The 
uadtree data structure seems to be the natural choice for the 
arallel architecture that has been used. The proposed approach 
Tovides the robust control over the mesh density. The basic 
'arallel architecture used for implementation is a network of four 
'800 transputers connected in a star network. The present 
rnp lementat ion generates a fully automatic mesh for two 
limensional geometries of any shape, singly or multiply connected, 
"he work done is in PARjC and provides some eKample results also. 
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CHAPTER 1 


INTRODUCTION 


1.1 Background Information 

1.2 Objectives, Scope and Limitation 

1.3 Organisation. 


1.1 Background Information 

The discretisation of the given geometry is widely studied 
problem as it has got wide applications. These problems are 
fundamental to many areas such as solid modelinSr computer aided; 
design, graphical rendering and scientific computations. 

A partitioning problem of particular interest, is optimal i; 
tr iangualt ion of a planar point sets. This problem finds ; 
applications in image processing, contour mapping, cartography, j 
finite element analysis. ■ | 

ii 

Finite element analysis is an important area of application, s 
The evolution of finite element method has been one of the most ; 
notable advances in engineering practices over last 25 years. 
Typical applications are in Aerospace, Civil and structural Engg, | 
steel and concrete bridges, soil mechanics and fluid dynamics. | 
Finite element analysis is com-monly used to obtain information i 
about structural behaviour without building prototypes. Today it 
is primarily used as an after the fact verification tool. Using 
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FEM in the earlier stages of the design can improve the 
design concept. Actual modelling and prototyping is not 
necessary. 

The following steps are followed for analysing the structure:- 

1. A number of attributes such as material properties to define 

physical behaviour under elastic deformat ion r loadcase 

includes applied forces and restraints and thus simulates a 
test case and future position of the part within the 
assembly. 

2. A mesh is generated from the boundary description of the 

part. This divides the interior of the objects into simple 

elements of the type known to the application program. In 
practice quadrilaterals or triangles for 2D and bricks, 

wedges and tetrahedrons for 3D are used. The mssh is 

a discrete approximat ion of the continuous interior of the i 
object. It is important that the difference between mesh ; 
and object should be small. The quality of the ■ 
appox imat ion is resolution. 

3. The loadcase information is assigned to elements and nodes [ 
and the mesh is analysed. The comiputat ions include matrix [ 

I 

approximations such as inversion or Guass type elimination. 
The cost of analysis is largely dependent upon the size of 1 
stiffness matrix which itself is a function of number of | 

f; 

degrees of freedom to which the accuracy is inversely | 

proportional and dependent further on mesh properties. 

The analysis is the discrete approximation of the 

continuous deformation of the part. Therefore more or 
higher order elements are required to model accurately 
regions of dramatic changes in the strain energy. Element 
numbering, global node numbering and nodal coordinate 



information together with additional element and nodal 
information such as boundary node data are required for 
these calculations. 

Thus inherent in all FEM methods is the necessity to divide 
the region of interest in a finite number of points or nodes. 
This subdivision is termed as mesh generation. This process plays 
an important role in the analysis of a design and also one of the 
fTiore sensitive stages in the process as it is intimately linked 
with both the geometry of the path and particular load case 
concerned. The accuracy of finte element solution is dependent on 
the element subdivision employed and generally makes largest 
contribution to the error involved. 

The basic criterion satisfied by any FE mesh subdivision is 
that piecewise variation of the basic unknown from element to 
element closely approximates the true solution. Finite element 
solution performed with increasingly finer mesh division 

essentially converge into the exact solution of the problem. But | 
it can be computationally expensive in the sense that larger size i 
matrix handling is involved with more number of elements- Thus I 
good balance of accuracy and cost is achieved. \ 

I 

The mesh generation is very time consuming and tedious jobrbut I 
is the essential part of the process. It can be error prone if 
done manual ly. This problem was identified long back and the need \ 
of automation was felt almost at the same time as of the i 
automation of FE method itself. Automation of the process j 

eventually relieves the user of tedious task of geometry 

definition and visualisation of the meshr which leaves him with 
the essential work for preparing the raodelrits attributes and 
specific loadcases. I 



Because of the ad vancement of the various CAD tools, it is 
possible to support visualisation and verification of the 
structure of nodes within the elements and of elements within the 
mesh . 


1.2 Objectives and Scope 


Faster processors and parallel architectures are becoming 
widely available. This development has motivated the present work 
to exploit the possibilities of using such architectures for 
faster FE methods. Though there are serious efforts towards the 
parallelisation of the FE method itself, as its inherent nature 
makes it attractive option far paral lei isat ion, there are very few 
attempts to parallelise the rriesh generation. So this work is the 
attempt to check such possibility. The advantages of parallel 
implementation are quite obvious. 

The objective of this work is to generate fully automatic 
mesh generation parallely. It is implemented on the parallel 
architectural class of tro.nsp-nt&z's. The approach used here for 
the mesh generation is the modified oxiadtree approach. 

There are various aprroaches, which are discussed in the next 
chapter, to the mesh genertion . But the quadtree method is an 

) 

attractive alternative in view of the paral lei isat i on of the | 
process and seemed to be the natural choice for the architecture I 
we used. The architecture we have used is a 4-node T800 tranputer 
network connected in star, as there is the minimum corrimunicat ion [ 
required arriongst the worker nodes. 

The system designed gives the' valid triangular meshes for mast of the: 



arbritarily shaped two dimensional geometries which are singly or 
multipaly connected. 

Although it has been limited to 2D domains it can easily be 
extended for 3D domains. Our aim is to provide a system designed 
to give the fully automated, faster, parallel model for mesh 
generation which gives the valid triangular meshes for most of the 
arbritarily shaped two dimensional geometries singly or multipaly 
connected. The work is limited rather to rriesh generation than to 
the actual Finite Element Method r equirenrient s . There are efforts 
to give the valid meshes as far as possible. But truly speaking 
the fully automted mesh gneration is not viable because in most of 
the cases the mesh generation depends upon not only the geometry 
but there are various factors such as boundary conditions, load 
distribution etc which really influence the final mesh generation 
requir ements. So features for interactvity should be provided for 
any practical mesh generators. 

1,3 Organisation 

The remainder of the thesis is organised as follows : 

In chapter 2 various mesh generation techniques, serial 
as well as parallel, have been discussed and corripared. It also 
contains the detailed- discussion about the modified quadtree 
technique which has been primarily used. 

In chapter 3, spatial hierarchial structures like 
quadtrees and octrees have been discussed. Various algorithms, 
serial as well as parallel have also been discussed. Its merits 
and demerits have been checked. 

In the last section, the actual implementation details has 


Deen provided. 



CHAPTER 2 


MESH GENERATION 


2.1 Introduction 

2.2 Serial 2D Mesh Generation 

2.3 3D Mesh Generation 

2.£. Pa>'sllel Mesh Generators, 

2.5 Modified Quadtree Approach 

2.6 Mesh Smoothing 


2. 1 INTRODUCTION 

Fundamental to the automation of FEM is the ability to 
generate a computat ional grid. It is the process of generating 
finite element models for simulated structural analysis. 

A fully automatic mesh generator is defined as the one which 
takes as an input ^a geometric representation and associated mesh 
control infartTiat ion and automatically generates the valid mesh 
without any manual intervention. A valid mesh is topologically 
cervpatible and geometrically similar to the model. Topological 
compat ibi 1 ity is a requirement that the mesh topology and 
geometric topology are congruent , while geometric similarity is a 
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requirement that the mesh geometry is an acceptable refinement 
which matches the geometric model exactly. Automation of the 
process relieves the user of tedious task of geometry definition 
and visualization of the mesh , which leaves him or her with the 
essential task for preparing the model rite attributes and the 
specific load cases. Mesh generation plays an important role in 
the analysis of a design and also one of the more sensitive stages 
in the process as it is intimately linked with both the geometry 
of the part and the particular load case concerned. 

The accuracy of the FE solution depends on the element 
mesh layout and the cost of analysis becomes prohibitively 
expensive if the number of elements in the mesh are too high. 

A good mesh generating scheme, therefore, should help the user in 
generating a mesh that is just fine enough for an adequate 
solution accuracy and efficiency. The conformity of the elements 
should also be assured to avoid gaps in the structure. 

The requirements of a good element are: 

1. Element should not be too elongated 

i 

2 Areas of steep gradient of stress are expected to have 
finer elements subdivision 

3. Quadrilateral elements do not offer much flexibility 
and mesh grading is a problem in these cases. So 
generally the triangular elements are preferred. Also 
in case of triangular elements compatibility is assured. 

4. The triangular elements should not be elongated wherever 
possible and triangles sho;jld not be obtuse 
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angle triangles as this produces stiffness matrix which 
may become singular or unstable 


Approaches to automatic mesh generation are many and 
varied substantially and are influenced heavily by the 
particular geometric modeler. The focus of mesh generation has 
moved from co-ordinate transformation scheme which is input 
sensitive and automatic tr iangulat ion which is fairly robust but 
involves time consuming searching to the hierarchical spatial 
decomiposit ion schemes which are versatile robust and are based on 
certain principles of computational geometry which enable them to 
function without applying too many heur istics. CCHE891 


2.2 SERIAL 2D MESH GENERATION 

There are serious efforts of automating the mesh generations 
part. There are over 250 methods and most of which deal with 
manually partitioning the region with quadrilaterals and 
occasionally to the triangular zones. Although there are lot of 
different methods adapted for serial mesh generation very little 
literature is available on parallel mesh generation. Whatever 
material is available on parallel mesh generators is from the view 
of parallel triangulation of given points? the technique which can 
possibly be used for mesh generat ion. CCHES93 
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Figure Classification scheme for mesh generation methods 


^ 'Table 1 . Comparison of the mesh generation approaclies 
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Current serial mesh generating schemes can be classified as: 

1. Interpolation mesh generation 

2. Automatic triangulation 

3. Quadtree/Octree approach 

4. Mesh generation based on Constructive Solid Geometry 

The interpolation mesh generation approach CCAV743par tit ions 
the structure into simpler subregions first and then meshes each 
subregion individually. Though this technique is quite popular 
because of its ability to produce well conditioned meshes. It is 
relatively difficult to carry out local mesh refinements because 
of the need to fit together in a consistent manner, the individual 
meshes defined by each subregion of the structure. 

Automatic tr iangulat ion [316783 has been popular method used 
in 2D mesh generation. However, this technique occasionally 
produces ill conditioned meshes and also it requires the eKtensive 
checking to ensure that newly formed elements do not pierce 
through or intersect any already defined ones. 

The third approach is based on Quadtree encoding of the 
object in 2D and Octree encoding in 3D eKtension. CYSHS3] The 

object is enclosed in a square universe which is then recursively 
subdivided into quadrants. The quadrants that are inside the 
object are retained and contribute to the elements of the mesh, 
those outside are discarded and those intersecting with the 
boundary are further subdivided. The subdivision continues till a 
prescribed level of the subdivision is reached. The quadrants at 
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the last level of the subdivision and intersecting the boundary 
are cut with respect to boundary and selective actions are taken 
to ensure that the result is a valid set of triangular elements 
All the in quadrants are also split to ensure a valid mesh. 

Mesh generation based on CSS (Constructive Solid Geometry? 
tree, generates the meshes for the object which can be represented 
by a CSG tree. CLEESSl It uses CSG tree to construct a set of well 
distributed points within the space of CSG objects and then uses a 
decision making process to construct a mesh of a this set of 
points. Point membership c lassi f icat ion is extremely useful in 
point generation algorithms. The efficiency of mesh construction 
depends upon the size of the decision tree. This approach can 
generate a good mesh for a CSG object with constant mesh density 
th^'oughcut the object. However, extension of this approach to 
variable mesh density needs special care. 

2.3 3D MESH GENERATION 

Standard element type in 3D Finite element method includes 
hexagon (brick), pentahedroni^edge') and tetrahedron. These 
building blocks must be assembled in a consistent way , i.e. 
vertex to vertex fashion. It is extremely difficult to decompose a 
solid into valid i.e. disjoint, compatible and consistent union of 
solid finite element elements. The rr*esh generation problem 
becomes especially difficult when a variation of element density 
from region to region is required. Local mesh refinement is often 
a rule than the exception. [CAV851 
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e lemertt 


mesh generat ion 


The three dimensinal finite 
techniques can also be classified as 

1. Interpolators 

2. Automatic Tr iangulat ion 

3. Semiautomatic Point Insertion Method 

4. Octree Modeling Technoique 

The most successful technique is called as isopararntr ic mesh 
generation. In this method an isoparametric matching is 
constructed fromi the conical domain (e.g. a unit cube) on a block 
structure. Then this mapping includces a natural coord i nat isat ion 
on the block with respect to which elemental node types can be 
referenced. This approach suffers from the serious problem of 
grading. CZIW833 

In the second approachr which is simmilar to 2D 
triangulation, the strategy of element generation is to select an 
edge connecting two node points in the solid and than surround 
this edge with tetrahedron defined by nearby points. During 
generation extensive checks must be made ’to ensure newly formed 
elements do not pierce through the early defined ones. Automatic 
tr iangulat ion is the same as that of 2D elements , the only 
difference is the elerrients are tetrahedron shaped . CCAVS5] 

The octree modeling technique is along the parallel lines to 
the modified quadtree technique. 

2. 4 Parallel Mesh Generators 

/ 

Although there is little disagreement for the advantages of the 
parallel methods few points should be considered before 
parallelisation of any method. Speed is the rriain advantage of the 



parallel architectures, but all algorithms can not be efficiently 
paral lei ised. Parallelisation is essential for the speed and it is 
important to explore the new areas of applications where parallel 
machines can be used. The method used here is an attempt to 
explore such a possibility. 


E. Mar ksCMARSSl gives the parallel algorithm for the triangulation 
of the given points. He gives the optimal parallel algorithm for 
triangulating a set of points in the plane. Given n points in the 
plane, join them by non intersecting line segments such that every 
region interior to convex hull is a triangle. The algorithm is of 
odog n> time using o<r)) processors and oCn) space. 


F. Chang et al.LCHES91, have proposed a new scheme which uses 
a vertex label assignments scheme to provide the information for 
the image generation so that the parallel ttissh generation on the 
basis o-*^ individual faces becomes possible. No checking of 
conformity is required as it is automatically assured. It carries 
the process of mesh generation by controlling the labels assigned 
to vertices rather than subdivision levels assigned to faces of 
regular quadr i lateral mesh. The subdivision process, which can be 
of two types, that is balanced or unbalanced is entirely driven by 
the values of labels assigned to the vertices of the input network 
the algorithm performs subdivision of all the faces 
simultaneously , each by one processor. Each face of the network 
represents the root of the quadtree and is subdivided and labels 
are also assigned to the vertices of these subquads. This process 
is continued till all the labels are equal to zero. 
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The parallel modified approach , which we have used is very 

attractive alternative for the parallelisation of the mesh 

generation method. We can see that very few algorithms are 
inherent nature of parallelisation. Quadtree is one such data 
structure which can be used for such purpose. Generally the 

communication overhead amongst the nodes decries the advantages of 
pa^'al lei isat ion. But in our case as each quad can be handled 
independent lyr minimum communication is required in between the 
nodes or server and worker model. In this method , each quadrant 
is handled independently by each processor .Col lecting the node 
information is only done by the server. The main drawback of the 
parallel programs is that they are highly architecture dependant. 
The parallel msodel of transputers is the attractive alternative 
foT- the quadtree generation which is explained in chapter 4 in 

details. The basic modified quadtree technique is explained below. 


2.5 MODIFIED QUADTREE METHOD OF MESH GENERATION 

The Modified Quadtree rriethod uses, conventional quadtree 
encoding technique with some modification for the 
tr ianguiar isat ion technique. CYSHS33 Mesh grading is proper in 
this case, in such a way that, a finer mesh in the regions of 
expected steeper gradient of field variables and a coarser mesh in 
the other places is possible. 

The Quadtree approach seems attractive because of its inherent 
nature of paral lei isat ion. As faster computers are available, so 
autoiTiating the preprocessor for FEM which is virtually the whole 
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time consumer and a tedious task, this method seems interesting. 

The Quadtree mesh generation schemes goes this way: 

In conventional coding of the quadtree, the object of 
interest is placed inside the square, that already has an integer 
coordinate system defined within it to identify the squares used 
in the quadtree representat ion. This square is subdivided into 
four quadrants. Next each quadrant is tested to see if it is 
inside the oh ject ( 5iL.4C/0 , or completely outside the object (.MftlTE} 
or only partially in . Homogeneous quadrants can be flagged 

as black or white , but the gray quadrants are again subdivided 
into four subquadr ants. These quadrants are tested in the same 
manner as above until some prescribed resolution level is reached. 
The advantages of the object representat ion by this way are, the 
storage requirement is far less than the pixel repr esentat ion and 
also the boolean operations, translat ional , rotational operations 
are faster. 

But the object represented in this traditional way has the 
drawbacks also. They are : 

1. In some cases the interior of the object may be represented 
by very large elements than the boundary elements- 

2. SofTiG neighbouring quadrants may undergo more subdivision 
than the others. 

3. Since the object is represented as the set of squares , non 
vertical and non horizontal boundaries had to be expressed 
in terms of the squares. So there is a chance of loosing 
some of the information while coding in this way. 



While encoding the object non horizontal and non vertical 
boundaries undergo repeated division engrossing the 
merory requirements. 

These problems can be solved by modifying the quadtree 
encoding technique. Quadtrees when modified use the less space 
which is the main objective. In this technique, there is one more 
type of quadrant defined as the cat qusdrant . This quad can have 
co'^'ners cut in any fashion. Thus when the required resolution is 
achieved, instead of representing it by hlsck nodes, it is stored 
as a cat quadrant, with the intersection points inf ormat ior. . In 

this way, we can represent the cut quad for any angular cut and 
there is no need of the extensive checking and the r ead j ustments 
G-f the cut points to boundary points as proposed by Yerry and 
Shepha-'d-'s scheme . CVSHeSl 

Ir. this way we ca^’ not have entirely integer tree , and hence 
there are no problems such as storage of the array and retrieving 
the information. Because in any case the floating point 
information has to be stored. 

The key to information stored at the node level are the 
discrete intersection of entities defining the boundary of the 
object with the boundary features of the node. The mesh is 

generated leaf by leaf method. The tree structure is used to 
transfer cell face by triangular information to neighbours to 
ensure the matching of the meshes within each leaf. The tree 
structure is also critical for containing the distribution of 
elements throughout the domain. The current implementat ion of the 



procedurE rsquires the neighbouring leaf nodes differ in size by 
only one level. This ensures the generation of smooth mesh 
gradation. Since the element rize is dictated by leaf size, any 
spatially based procedures that can dictate the size of cells, can 
be used to control the distribution of the elements throughout the 
mesh. 


Use of an optimal mesH should produce a 
satisfactory accuracy for less computational cost, 
optimal tneshes are extremely difficult to construct 
the cases. So in this case we can start with 
solution of coarser grid and develop procedures for 
iiTiproving it by introducing finer meshes in the 
accuracy is inadequate. 


solution of 
Unfortunately 
for most of 
a preliminary 
automat ical ly 
region where 


Just construction of Modified Quadtree do not produce 
satisfactory element mesh since the quadtree/octree approach 
itself is not an element generation technique, rather it serves to 
simplify the geometry for an associated element generation 
scheme. The rich data structure provides a high degree of local 
control over the meshing process. So additional features must be 
added to create transitions between quadrants of different levele 
to ensure reasonable aspect ratio and obtain beet possible mesh. 
To produce a valid mesh in which no common nodes lie along with 
the edge of the another element, transitic" meshes must be 
replaced between the quad’^ants of different levels. But in order 
to limit the number of multilevel transiticm meshes, e second 
modification is carried out on the Riodified quadtree to ensure 
only one level difference between any two neighbouring quadrants. 



Afterwards the internal and boundary nodes can be triangulated 
independently to give valid mesh. The element are improved by the 
application of Laplscian operation. This operation places the 
interior nodes at the centroid of the nodes that they are 
connected . 


The finite quadtree and octree mesh generation procedures are 
well suited for the adaptive mesh refinements also. Starting with 
an initial coarse mesh the nodes which are having large gradients 
can be bisected independently and the new mesh is generated. 
Conputat ional growth has been used as a formal indicator of the 
performance of mesh generations schemes. THis indicator measures 
the rate of increase of CPU time with increasing number of 
elements. The advantage of this measure over the use of just CPU 
time alone is, it is independent of the computat ional speed of the 
machine used. In general a plot of CPU time vs number of elements 
is used as a measure c^ growth rate of mesh generation scheme. 
The best growth rate upto this time is found in quadtree schemes 
in which the CPU time is linearly proportional to the number of 
elements generated. 

quadtree mesh generation scheme localises elements formation to 
within the individual quadrants. When the information must be 
exchanged between a terminal boundary quadrant and its neighbours 
there are only four directions in which to conduct- a search, no 
matter the level of the quadrant. More over the interior 
quadrants can be directly triangulated with out any further 
processing. This results in a linear rate for the quadtree mesh 
generation scheme. 
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2.6 MESH SMOOTHING 


After the tr iangulat ion process the mesh smoothing techniques are 
applied. These are repetitive processes which reduce the overall 
speed of the execution. Often the elements produced by an 

automatic mesh generator are not well shaped. In such cases it is 

« 

necessary to apply the mesh smoothing techniques to improve the 
mesh. The most popular technique is Laplacian technique which 
seeks the repositioning of of the nodes such that each internal 
node is at the centroid of the polygon formed by the its connected 
neighbours. The repositioning is usually done iteratively. 

Automatic mesh generation has still number of disadvantages. 
No CAD system can supports necessary handles for the mesh 
generation which will allow material properties and load cases to 
be specified with the geometric model. All automatic mesh 
generators are therefore, forced to create more dense meshes than 
would be otherwise be necessary. In order to replace the mes*^ 
generation schemes by more reliable scientific methods and total 
automation, a more precise knowledge of convergence and absolute 
errors of results are necessary. 

The state of the art in mesh generation is such that fully 
automatic mesh generation with adaptive capabilities is not yet 
fully a reality. The requirement of mesh density is based on the 
observation that the use of less dense mesh either locally or 
globally after rediscretisation can lead to poor results with the 
same problem solved with remeshing. 
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CHAPTER 3 


QUADTREES AND OCTREES 


3- i Introduction to quadtrees and octrees. 

3.2 Parallel irriplementation o-P quadtrees. 

3.3 Octrees. 

3.4 Algorithms for Quadtrees. 


3. 1 Quadtrees and Octrees 

There are various data structures to represent image data. 
There are numerous hierarchical data structuring techniques in use 
fo»' representing spatial data. Hierarchical data structures a^'e 
important techniques in the areas of image processingj- computer 
graphics, solid modeling, geographic information systems etc.. 

They are based on the principle of recursive subdivision. They 
are mainly useful because of their ability to focus on the 
irfteresting subset of the data. The hierarchical organization of 
the image space allows the resolution to vary with the complsKity 
of the objects in various regions. They are used as a basic 
geometric modeling data structures. Many operations on the 
images can be performed with other types of data structu’'es, 
however, hierarchical data structures are attractive because of 
their conceptual clarity, ease of implementation and inherent 



character ist ic for paral lal izat ion. There are various types of 
hierarchical data structures as pyramids, quadtrees, octrees, 
polytrees. The term Quadtree is used to describe a class of 
hierarchical data structures whose common property is that they 
are based on the recursive subdivision of the continuum. It is a 
type of data structure used for compact representation of 2D 
images. It is based on the recursive subdivision of the 

space. Quadtrees are same as the normal binary trees but with 
degree four. Generally a node is either a leaf or having the four 
sons. A quadtree is generated by dividing the image with 
quadrants and repeatedly subdividing the quadrants into 
sub_quadrants until each quadrant has uniform color. The stopping 
rule of the formation of the quadtree is a leaf criteria- The 
iTiodi f icat ion o^ leaf criteria can be the depth of the tree. Once 
the may.irrium depth is reached , even if the leaf criteria is not 
satisfied the region is simply represented as a leaf nods. 

But in such a case it is possible to loose some of the 

informiat ion. The leaf criteria may be that the quad may lie 

completely inside the object, represented as BLACK node or leaf, 
or it may be corripletely outside the object, represented as Mf/ITE 
node, or it may be cut at the boundaries of the csbject, 
represented as GRAY node, for which -again subdivision is 

necessary. Thus heterogeneous node has four descendants called as 
Sfe', S£, Nt/, N£ corresponding to the southwest, southeast, 
northwest and northeast quadrants of the parental quad Thus a 

node at level h corresponds to 2 sized quad, where 2 is the 

size of the original image space. This can be modified to the 
maxirrium depth technique. In this technique ,Hven if the leaf 

criteria is satisfied, subdivision takes place of completely 
inside nodes till the maxifrsum prescribed level of the tree 

formation is not reached. 



Figure . Pointer encoding of the Quadtree of Figure 
. . Internal nodes are represented by circular nodes. 
Terminal nodes are represented by square nodes 
whose contents correspond to the blocks in Figure 



The quadtrees can be di f ferent iated on the following basis 

1. The type of data it represents such as region, point, 
curve , surface or voluiTse. 

2. The principle of decomposi t ion process such as regular 
or irregular quadtrees. In case of regular quadtrees 
decomposi t ion at each level is into the equal parts 
such as hexagonal , square quadtrees tBAI1SS3 

Regular decomposition type of quadtrees are easy to formulate 
and program and not much heuristic is involved. In case of 
irregular division, the division is governed by the nature of the 
input data. Even though there can be compact representat ion, too i 
miuch heuristic is involved in this. There are various ways to 
store the quadtrees, such as classical quadtrees, linear 
quadtrees. The amount a* storage required by the quadtrees or 
octrees is directly proportional to the number of lea-? nodes. 
Generally a node is represented by a record with pointers to the 
four children. There is also a field for the identification of the 
type of node such as 5r-aV, black or iiihite. If the node is a leaf 
nods than all pointers to son aye NULL. The main limitatiori of 
quadtrees is that the the norrrial classical structure makes it 
spacHwisH expensive as the leaf nodes have the null children. So 
various other ways of storing the quadtrees has been explored. 
Linear quadtrees is one of the irjodi float ion for the repressntat ion 
of quadtrees. In this method, the pointer less repressntat ion of 
the quadtrees has been implemented. But classical quadtrees are 
conceptually clea'^ and traversal algorithms are very easy and 
faster. Mast of the quadtree algorithms are simply preorder 
traversal and hence their execution time is a linear function of 
number of nodes in the quadtree. The root of the quadtree 
corresponds to the image it represents. A node in a quadtree is 



either a leaf or has four node son i.e. non terminal node. Each 
son node is associated with a quadrant of the block corresponding 
to the father node and the original image. CSAMS41 ** 

The advantage of the quadtree representation for images is that 
simple and well developed tree traversal algorithms allow faster 
execution of certain operations such as superposition of image 
area and perimeter calculation. In addition the 2D coordinate of 
each block is implicitly stored in and can be readily recovered 
from the quadtree representation. So there is no need of 
explicitly storing the information of the co-ordinates of the 
points of the quad for which the node is represented. Each node of 
a quadtree represents a definite corresponding block in the 
original image. Given a polygon with p vertices , construction 
of quadtree within the space of 2 k 2 has execution time of 
0(v^n>. One of the most useful features of the quadtree is that 
an increase of the resolution parameter requires further 
subdivision of the terminal nodes only, leaving the remaining part 
of the quadtree unchanged. Roughly speaking a refinement of the 
image merely requires the refinement of the leaf only. Thus 
quadtrees are dynamic structures. 

3.2 Parallel implementation of the quadtrees 

There are various parallel quadtree building algorithms. Dyer 
has described a bottom up algorithm which builds a quadtree by 
recursively merging the blocks of the same color starting from the 
pixel level . tDYRSll. Eiedeman has presented a Concurrent Prolog 
algorithm for quadtree building. EEDRS51 Another approach is using 
shuffled row major indexing scheme for the pixels in the given 
image. Here, for each pixel it is determinedr whether it cars 



represent a leaf node and if ,it can ,the size of the leaf node is 
determined by the number of its preceding pixels and the number of 
its succeeding pixels in the run it belongs to. The trend in the 
parallel algorithms to construct the quadtrees is mostly from the 
image processing point of view. These algorithms are basically 
meant for the data which generally comes as a raster or pixel 
format. They are mostly bottom up approaches in which the tree is 
formed by merging the same coloured neighbouring pixels to form 
the leaf. The approach generally changes depending upon the format 
in which data comes and how it is arranged. There are not many 
explorations such as the top to bottom approach and any other 
ways. CS83 

3.3 Octrees 

The modeling of 3D objects and space by computer is irriportant 
in space planning , computer animation and a machine vision. An 
extension of the quadtree structure is tfie octree structure for 
the representation of the 3D objects. The octree generation is 
done or. the sartie basis as the quadtrees. An octree is generated by 
dividing a universe into octants and sub-dividing the octants into 
sub-octants with all voxels in each octant lie entirely within or 
outside of an object. Thus octree nodes which are called as voxels 
represent the cubic volume in the original image. The voxel is 
either a leaf or internal voxel. Internal voxels have the eight 
children and ths Te-af nodes do not have children at all. 

The octree advantages over other methods for solid modeling 
are; Any arbitrarily shaped objects convex. concavEr with 
interior holes can be represented in the required precision. Only 



Figure . (a) Example 3D object, (b) its octree block 
decomposition, and (c) its tree representation. 








a single set of manipulation and analysis algorithms are required 
for all objects. New techniques are not needed to handle 
sophisticated or complex shapes- Because of spatial sorting and 
uniformity, hidden surface display algorithms requires no 
searching and sorting. The main advantage of the octree model is 
the linear complexity of the boolean operations and rendering 
algorithms. Several geometric properties such as surface area, 
interference are easily calculated- However they are approximate 
models as the slanted surfaces must be represented by the minimum 
scale black nodes. The modified quadtr ee/octree approach reduces 
the degree of sub_division and therefore requires less storage 
than the classical octrees. 

The main disadvantage is the large memory r equiremient for 
tncodihs. Th® algorithms generally are 0( surf ace area of the 
object) for the volutTi-etric data and for quadtrees generally 
0 ( perimetere of the object) 

There is a large literature available on various operations on 
the quadtrees such as 

1. Algorithms for the quadtree formation. 

2. various operations such as transf ormat ion or rota'tion of 
the quadtree. 

3. Neighbour finding. 

4. From raster image point of view. 

Out of that from this work point of view very few algorithms 
are useful. The one technique we have used is explained in detail 
over here. 



3.4 Neighbour finding algorithms 


In several applications, the computation of each node of the 
quadtree is dependant on the value of its adjacent neighbours. 
Thus we must be able to locate the neighbour of a node. 

Two tiles are neighbours if they are adjacent either along an 
edge or vertex. The requirement for our algorithm was to find the 
neighbour along the edge and not along the corners. The neighbour 
finding operation must performed in a most general way i.e. the 
neighbours must be located in a way that is independent of both 
size and the position of the node. Locating adjascent neighbours 
in the horizontal or vertical direction is relatively 
straight f orward. The basic idea is to ascend the quadtree until a 
common ancestor with the neighbour is located. And then descend 
back down the quadtree in search of the neighbouring 
node. Neighbours along the edges need not correspond to the block 
of the same size. If the neighbour is larger then only part of the 
path from the nearest comimon ancestor is retraced, otherwise the 
neighbour is of the same size of the node. If there is no 
neighbour as in the case of border nodes then NULL is returned. 
The nearest ancestor which is comman is the first ancestor node 
which is reached via its sons. We retrace the path used to locate 
the nearest common ancestor except that we make a mirror itriage 
moves about an axis fortried by the common boundary between the 
nodes. For example in the case of W neighbo*' the mirror ifnages of 
NW and SW are NE and SE. 



CHAPTER 4 


SOFTWARE DEVELOPMENT 


4. 1 Introduction 

4.2 Quadtree Generation 

4.3 Validation o" The Quadtree 

4.4 Node Numbering And Connectivity 

4.5 Parallel implementation 


4. 1 INTRODUCTION 

Here in this chapter , actual implementat ion details are given 
Also given is how the modified quadtree method used for this 
purpose is different from the original method proposed by Yerry 
and Shephard's idea of integer tree. The details of the parallel 
architecture used are also provided. The parallel algorithm 
for mesh generationr implementation problemsr the lirr.itation of 
the architecture for the method to parallelise are also discussed. 
The next chapter gives the detail information of the results^ 
Comparison between serial and parallel methods and the scope for 
the future work. 

The actual implementation of the algorithm can be divided into 
five modules 



1. Data Acquisition. 

2. Quadtree Generation 

3. Valid Quadtree Formation. . 

4. Actual Tr iangulat ion. 

5. Global Node Numbering and Connectivity . 

The only input required for the modified quadtree technique is 
the boundary points. The method assures the valid mesh in most' of 
the cases, as the finer elements in the area of steeper gradient 
of variables and coarser elements elsewhere. So the minimum mesh 
refinement is required. But sometimes it is difficult to shift 
the newly created points on the boundary after foriiiation of the 
mesh. The drawback is particularly true for current method 
because of its spatial sensitivity. So sometimes required number 
of points on the boundary are also specified in order to get more 
refined mesh. Although we have assumed only straight line 
polygons as the input for the present work, it is easy to write a 
preprocessor which can accept any shaped geometrical object 

input, from any standard CAD package and is suitably converted 
into its polygonal representat ion. 


4.2 Quadtree Generation 


Once the input polygonal data, in which polygons can have any 
shape, either singly or multiply connected, is acquired, it is 

n rt 

enclosed by the square o-^ size 2 » 2 n depends upon the 

co-ordinate information of the polygon. As already explained the 
quadtree encoding goes like this as : 



The original square is divided into four equal subquadrants. 
Now each quadrant is checked against polygon whether it 
completely inside the object (BLAC/O, or completely outside the 
object or partially cut by the polygon boundaries. 

The GRf^y quadrants are again subdivided in the same way. The 
quadtree formation is controlled by two parameters , M3x_depth 
and Min-ResolutioT), Max_depth gives the maximum depth of the tree 
upto w’^ich the tree formation should gor irrespective of its 
status. Min_rHsolut ion defines the stopping criteria for the 
leaf. In generating quadtree even if the quad is 
completely inside the object , it is subdivided into four BLACK 
nodes if the required depth is not yet reached. The gray quadrants 
are termed as graV or cat according to the depth of the tree 
that has been reached. Instead of passing the original polygonal 
structure every time to each quad, only the parental cut quad 
polygon information is passed to its children. As the 
determination of the status of the quad is a function of number of 
edges cf the of the polygon, significant speedup is achieved. 
Normal classical way cf pointer tree has been used to represent 
the quadtree structure mainly for two reasons. One, as the ease of 
itTiplementat ion and maintainance and also it is conceptually very 
simple. Second the efficient neighbour finding algorithm 
explicitly needs the quadtree structure in the classical tree 
form. With other types of representations explicit pointers to 
nighbouring nodes are required, which kills the very purpose of 
using the other type representations. The main constraint to the 
classical quadtree representation is the pointer overhead. The 
parallel transputer model which we have used limits no way this 
representat ion as the individual nodes have sufficient meiTiory. 
(AHb RAM for each node). 



The individual status of the quad is determined against esch edge 
of the polygon. We have used Meil er_fithertror> algorithm to 
determine the cut surface of the polygon against the 

quadrant. CWE1771 Although it is somewhat complicated to implementr 
the popular Sutherland Hodgensn algorithm is not used as it gives 
a lot of dangling regenerative edge cases for the polygons with 

the holes which are difficult to eliminate. CSUT741 

The main difference between Yerry -Shephards modifierd quadtree 
algorithn-i and the present work is that there is no concept of 

integer tree involved. The integer tree is unnecessary in the 
sense that final point representation should be, anyhow in terms 
of the floating points. So all basic purpose of using integer 
tree is lost. In the above method cut quadrants are represented 
as the integer points, in which the quad is subdivided into 
definite intervals. To determine the nature of the cut quad, 

almost eighty enumerated cases CGANS82 have to be explicitly 
checked. Also we have to store this information in a staclr which 
requires another complicated algorithm to retrieve the 

informat ion. Due to this also, after forming the mesh, there is a 
requirement to pull back the node points to the boundary anc! to be 
stored as a proper floating point representation. All "this can be 
avoided with a simple memory overhead of one pointer for each 
leaf, which is used to represent the resultant cut quad pol-y after 
cutting the polygon by the quadrant. There is no need to check 
the cases, retrieving the information and pulling back the 
boundary nodes. Also there is no explicit need to determine the 
boundary nodes as a complex boundary discrimination algori-thm is 
used in the above mentioned case. 



4.3 Validation of the Quadtree 


Once the quadtree is formed ^ it can not be used directly for 
the mesh generation. Additional features must be added and checked 
to create transition between quadrants of different levels. To 
produce a valid element mesh in which no corner nodes lie along 
the edge of another element. transition meshes roust be placed 
between quadrants of different levels. Such a situation allows 
gaps to open in the structure, unless the continuity of the 
displacement is ensured. So in order to limit the number of 
multilevel transition meshes, a modification of tree structure is 
ca^'ried out which ensures only a one level difference between any 
two neighbouring quadtrants in all four directions. In this case, 
the conformity of mesh should also be assured to avoid gaps in 
the structure, i.e. the intersection of two non disjoint, non 
identical elements consists of a common vertex or edge. The 
algorithm explained in the last section is used to find out the 
neighbours of each leaf in the four directions. If the level 
difference is found to be more than 1 then the particular leaf is 
individually subdivided without modifying the rest of the quadtree 
structure. Thus after checking the valid level difference between 
each node and all its neighbours as less than or equal to one, 
the tree is ready for direct triangular isat ion. 


4.4 Tr iangular isat ion. Node Numbering and Connectivity 

In this work only triangular elements are implemented since 
these are the most flexible elements. The individual leaf nodes 
tr iangular isat ion is carried according to the geometry it 
represents. The square nodes are represented by the diagonal 
elements, triangular nodes are accepted without any modification. 



Only the nodes with cut status and the nodes whose neighbour is at 
one level lower has to be carefully triangulated. 

Once the tr iangular isat ion pattern is decided , global 

numbering is given to the each point and global connectivity 

information is stared. Finite element also requires the 

connectivity numbering of the nodes in the specific fashion so 

that the connecting points should have minimum difference between 

their numberings. How to do this is itself an independent 

problem, for which, in this work has not been given enough 

attention. But a postprocessor can be written to ensure the 

minimum width stiff matrices. The elements aspect ratio can 

not be ensured in the elements produced in this way so these are 

irrsproved by the application of Laplacian operation. In this 

operation the position vectors Pi of the interior nodes are 

positioned to satisfy the expression P. = l/n. « EP- i=lr m 

^ 

where n. is the number of nodes connected to the node i, P. is 
1 J 

the position vector of the two connected nodes and m is the total 
number of node points. This equation places the node at the 
centroid the the nodes that they are connected to. 


4.5 Parallel implementation of Algorithm 


Although the advancement of technology has made the parallel 
architectures widely available, there are certain serious 
drawbacks generally faced by the parallel programmers. 
Existing parallel languages are tied to some architectural 
class. Because of paradigms and pattern of various parallel 
architectures differ, the programmer roust be aware of the kind of 
architecture on which the software is going to be run, the 
number of processors, storage capacity and their configurat ion. So 
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before representing the actual parallel algorithm the parallel 
architecture is briefly explained here and afterwards the actual 
parallel algorithm is provided. 

Parallel archi tectures can be broadly classified as 

1. shared memory models 

2. Message passing models. 

In case of shared memory models the comtrton memory is shared 
between all of its nodes and the issues like exclusive read/write 
operation comes into picture. In the other model of message 
passing there is no common memory between individual nodes. Each 
node has its own separate memory area, the communication is 
strictly through the message passing over the channels. Each model 
has their merits and demerits. The parallel architecture, 
which has been used is, transputer based. It is based on the 
message passing architectural model. There are four fSSOO nodes 
connected in a star fashion. Each node is having its own memory of 
4Mb RAM. There is a one node which is connected to host PC and is 
a master to which all the other nodes of the network are 
connected. The work done is in P^R_C. 

The treatment of parallel processing in transputer systems is 
based on the idea of cammunicating seQaential Processes. CH0R781 
In this model, a computing systerri is a collection of concurrently 
active sequential processes which can only communicate with each 
other over channel. A channel is nothing but the one sided 
communicat ion way which connects exactly one process to another 
process. Each process can have any number of output or input 
channels but they can not be dynamically created. The channel 
communicat ion is synchronised. A process wanting to send the 
message over a channel is always forced to wait until the 



receiving process reads the message. A process in this model is 
just a black box connected to the outside world only by its 
channel Sr thus the actual nature or implementation of the 
process is not important. Each transputer processor has four 
Inaos links to connect it with other transputers. Each link has 
two channels, one in each direction. These channels behave exactly 
like software channels and provide synchronised , unidirectional 
communication. Each processor can not be connected more than the 
four other processors. It is also possible to run any number of 
concurrent processes to be run on individual processor. The 
equivalence of internal channels means it is possible to develop 
a complete parallel system on a single transputer and then move 
some of its processes onto another processors without having to 
recompile any code. 


Parallel C is based on the same above described model of trans- 
puters. A complete application is viewed as a collection of one or 
more concurrently executing tasks. Each task has its own memory 
for code and data, a vector of input ports and output ports. The 
code of a task is a single transputer image file generated by a 
linker. Each element in a input and output vectors is of type 
pointer to channel word CCHRN *} Ports are bound- to- r-eal channel 
addresses by configuration software external to the task. This 
software also provides ways of specifying which software tasks are 
to be run on which hardware processors. 

For many parallel computations it is useful to be able to 
create applications which will automatically configure themselves 
to run on any network of transputers. Such applications will run 
automat ical ly faster when more transputers are added to the 
network without reconf iguration or recompi lat ion . This technigue 





is called as Processor farming. In the processor farm technique, 
an application is coded as one master task which breaks the job 
down to small independent work packets which are processed by any 
number of anonymous worker tasks, All worker tasks must run the 
same code. Each worker simply accepts the packet over the network 
processes it and sends the result packet back to the master via 
routing software. 

The quadtree method is attractive from parallelisation point of 
view. The architecture we have used is suitable for this. The 
model consists of a root node which is connected to the host of 
the PC. The other three node transputers are directly connected to 
the root. The communication between individual nodes is 
unnecessary for this application, since most of the information 
processing is done individually. Benerally speedup achieved by the 
parallelisation is marred by the heavy communicat ion requirements 
between the neighbouring nodes and the expected speedup is not 
obtained, but in this case there is a little communication 
overhead. In our case the communication overhead was minimum 
since we can build the quadtree truly independent of each quads. 
Only communicat ion to master is required. 

The five modules which should be carried sequentially, the 
paral lei isat ion is done according to the individual model. 

1. Data acquisition is done by the root node. 

2. The two tasks are there. One task which determines the 
status of the quad and the resultant cut polygon for the 
node. Other task is a master task which waits for the 
four workers to complete the work and takes the 
decision of the further subdivision of the quad. The worker 
task is kept on the four nodes. On the root the two task 



run concurrently 


Thus serial tree formation is a function of total number of 
nodes in the tree , while parallel version is a fttotal no of 
GRAY nodes in the tree) so for the larger depth of the trees the 
speedup of almost the number of processors is achieved. Thus if 
N is the total number of nodes in the tree. The number of Brav 
nodes = (N-i)/4 . Since the quadtree node has either four children 
or none. Thus a speedup of the factor 4 is achieved. 

The formation of the quadtree thus can be done parallelly. 
the validity of the tree has to be checked by the root node as it 
involves the information of the neighbouring node, which can not 
be done truly independently , The level difference between the 
node and its neighbour is checked parallely and the worker sends 
this information to the root node , where root node takes the 
proper decision of balancing of the tree and ensures the level 
difference of one in between all the leaf nodes and their 
respective neighbours. 

The valid tree leaf nodes can individually triangulated 
according to the type of the polygon it cuts and this information 
is passed to the root node . Root node accepts the connectivity 
information and the co_ordinates of the nodes and gives the global 
numbering as well as stores the connectivity information for the 
mesh generator's final output. 

1. Parallel Quadtree Generation 
PARJBUILD_TREE (tree, quad) 
if (tree->coIor == "GRAY^) then 

PAR_DET_COLORJDF_CHILDREN(trHB, quad)? 
subdivide quads into four children quads; 

PAR BUILD TREE(tree->ch!:i], quadl) 



PAR_BUILD_TREE(tr ee->chC2] , quad2) 
PAR_BUILD_TREE(tree->chC33 , quad3) 
PAR_BUILD_TREE( tree->chC43 , quad4> 


The procedure PAR JDET_COLOR_OF_CHILDREN is a task on each worker 
processor, altogether the four worker processes determine the 
status of the four eon quads. 

2. Validity of the tree 
For leafjnodes 

do 

assign leaf to each worker. 

GET_P ARJ-EVEL_D I FF 1 1 ea f ) ; 
if (levei_diff > 1) then 

PAR_DET_COLOR(lBaf, quad!; 

3. Generation of the Nodal connectivity 

assign each worker the leaf; 
get the local nodal coordinates and the connectivity 
informal ion. 

give the global nodal numbering and the connectivity 
inf or mat ion. 

The smoothing operation is handled by the root serially as the 
shared memory is not available. When the co_ordinate information 
is changed , the new information must be used by each worker . The 
model of message passing does not allow to do it efficiently as 
large communication overhead is involved. 

The next section gives the actual experimental results for 
serial as well as parallel implementation . Also proves the 
results experimentally. 
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CHAPTER 5 


CONCLUSIONS 

The algorithm described in the preceeding chapters has been 
successfully implemented on TRS 800 transputer based PC svstep. 
It has been tested on several data sets. The algorithm gives 
satisfactory results for a number of geometries such as star 
shaped polygonsr concave as well as singly or multiply connected 
polygon boundaries. There is almost no user intervention needed. 
The algorithm gives robust control over mesh generation using two 
parafTiGters: ilthe Ma>;_depth ii) Min_resDlut ion . 

Despite the fact that the algorithm works with different kinds 
of data sets it suffe'^s certain I imitations. The quadtree 

structure is hig'-ly sensitive tc the spatial arrangements. It is 
possible that for the sanve geometric structure defined in 
different co-ordinate references one may getdi-^^-^Erent 

tr iangulations. This is not an algor ithttfic problem but 

implementat ional difficulty. The geometric problerrfS are 
difficult to imple.ment since the number of possible cases is 
inordinately large. very di-*^-'^icu!t to consider all the 

specialised cases. In formation of the quadtree if such unstte^'ded 
cases arrived ^ then it is possible to get the erroneous results. 
For the present work such cases are avoided. 
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FIG. 3 GEOMETRY WITH HOLES 
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